Qubits entanglement dynamics modified by an effective atomic environment 
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We study entanglement dynamics of a couple of two-level atoms resonantly interacting with a 
cavity mode and embedded in a dispersive atomic environment. We show that in the absence of 
the environment the entanglement reaches its maximum value when only one exitation is involved. 
Then, we find that the atomic environment modifies that entanglement dynamics and induces a 
typical collapse-revival structure even for an initial one photon Fock state of the field. 

I. INTRODUCTION 

Entanglement generation between two atomic qubits has attracted considerable attention during the last two decades 
due to its importance in various quantum information processes [lj, Q- Those ideal processes, such as quantum 
teleportation, quantum cryptography, and quantum computation algorithms are strongly related to the capability of 
generating bipartite entanglement However, in real quantum systems there are uncontrollable interactions 

with the surrounding environment which usually lead to a decoherence resulting in the destruction of the entan glem ent. 
Recently, several effects of different kinds of noisy environments, specifically bosonic environment |fi r7l l8t Ifl llfl lllj. 
and fermionic environment jl^. Ilil IhH. llfl Il9j| on the entanglement dynamics have been extensively studied. Especial 
effort was applied to find decoherence free entangled states [14j,ll5|- For instance, B. Kraus and J. I. Cirac $] show 
that two atoms can get entangled by interacting with a common source of squeezed light and the steady state is 
maximally entangled even though the modes are subjected to cavity losses. S. B. Zheng and G. C. Guo (2l| proposed 
a scheme to generate two-atom EPR states in such a way that the cavity is only virtually excited. S. Schneider and G. 
J. Milburn |ll| show how the steady state of a dissipative many-body system, driven far from equilibrium, may exhibit 
nonzero quantum entanglement. Molmer and Sorensen have proposed a scheme for the generation of multiparticle 
entangled states in ion traps without the control of the ion motion. 

Although the effect of the environment on the atomic entanglement is usually destructive, in some specific situations 
two quantum systems can get entangled in the process of their decaying to a common thermal bath [16L Il7| . A similar 
effect was discussed in 0| where a method of generation of entangled light from a noisy field has been proposed. It was 
also shown |19| that the interaction between two spins and an itinerant electron environment leads to entanglement 
of the initially unentangled spins. 

In this article we study how an effective atomic environment modifies the atomic entanglement generated in the 
course of resonant interaction of a single mode of the cavity field with a couple of two-level atoms (the so-called Dickc 
or Tavis-Cummings |2 ll model) . Evolution of entanglement in the two-atom Dicke model was previously studied in 
the case of an ideal cavity in |22j| and in the presence of a dissipative environment in . Our study is motivated by 
the following physical situation: consider a cluster of two-level atoms (resonant with a mode of a cavity field) placed in 
a strong electric field (see e.g. H3,|2^|). Physically it could be a cluster of polar moleculae. The electric field generates 
a noticeable Stark shift so that most of the atoms are detuned far from the resonance, except a very small portion of 
them, whose dipole moments are approximately orthogonal to the field. Because the atom- (quantum) field interaction 
times are much shorter than the typical times of atomic diffusion, we can consider that the orientation of the dipole 
moment is "frozen" and that the physical mechanism of changing the atomic dipole orientation is a collision with the 
cavity walls, since collisions between the atoms in an atomic cluster are practically improbable. In the process of 
interaction with the cavity field the resonant atoms become entangled. We will study the simplest situation where 
there are only two resonant atoms. Nevertheless, the effect of non-resonant atoms on the dynamics of resonant ones 
is not trivial. The dispersive interaction of the field mode with non-resonant atoms leads to a modification of the 
field's phase which, in turn, affects the evolution of resonant atoms. Thus, the non-resonant atoms play the role of an 
effective dispersive environment whose whole effect could be expected to reduce to a phase dumping |26j , and thus to 
the entanglement decaying. Nevertheless, as it will be shown, the influence of such effective environment is not always 
destructive but also leads to a constructive interference, which reflects in, appearance of a system of collapses and 
revivals of the atomic concurrence even in the presence of just a single photon in a cavity. The article is organized as 
follows: In section [H] we analytically show, for some specific initial conditions (non-excited atoms and the field in a 
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Fock state), that the entanglement of formation in a bipartite system of two-level atoms interacting with a quantized 
mode reaches its maximum value when only one excitation is involved and it decays as 1/n when n» 1, being n the 
number of photons in the initial Fock field state. In section ITTT1 we derive the effective Hamiltonian of noninteracting 
two-level (resonant) atoms and a cluster of A atoms (far from resonance) interacting simultaneously with a quantized 
mode and we find the evolution operator when only one excitation is considered. In section llVl we study the effect of 
the dispersive atomic environment on the entanglement dynamics generated by one excitation for two different initial 
conditions. In section we summarize our results. 



II. ENTANGLEMENT IN THE TWO-ATOMS DICKE MODEL 



By entanglement of two subsystems we mean the quantum mechanics feature whose state can not be written 
as a mixed sum of products of the states of each the the subsystems. In this case the entangled subsystems are 
no longer independent even if they are spatially far separated. A measure, E(\(f>}), of the degree of entanglement 
for a pure \<p) state of a bipartite system can be given by means of the entropy of von Neumann, of any of the 
two subsystems. For a mixed state p the entanglement of formation E(p) between two bidimensional systems is 
defined as the infimum of the average entanglement over all possible pure-state ensemble decompositions of p [30| . 
Wootters found an analytic solution to this minimization procedure in terms of the eigenvalues of the R = \J y/pp^fp 
or R! = pp non-Hermitian operators, where the tilde denotes the spin flip of the quantum state. The solution for 
the £(/o) concurrence associated with the entanglement of formation of a mixed state of a bipartite of bidimensional 
subsystems is given by €.(p) = max{0, Ai — A2 — A3 — A4}, where the A^'s are the square roots of eigenvalues of the R' 
operator and the eigenvalues of the R operator, decreasingly ordered. Throughout this article we consider this £(p) 
as a measure of the entanglement degree between the two resonance atoms a and b. 

Let us consider two identical two-level atoms resonantly interacting with a single-mode cavity field. The interaction 
Hamiltonian has the form 

H = g{a(s +a + s +b ) + a f (s_ a + s_&)}, (1) 

where s+j = \l)jj(0\ and S-j = |0)jj(l|, with and \0)j being the excited and ground eigenstates of a Z j of the 
jth atom (j = a,b), a and are, respectively, the creation and annihilation operators for the cavity mode, g is the 
atom-cavity coupling strength. Considering initially a Fock field state and both atoms in their ground states, the 
reduced atomic density operator, at time t, is 



[nC n (t) +n- l] 2 nSl(t) + , n(n - l)[l-C n (f)] 2 
Pah = x y |0) o |0)& o (0|&(0| + 2n _ l \V )ab I + (2n- l) 2 WaWbaWbiM, (2) 

where we have defined the functions: 



C„(t) = cos(v/2(2n- rjgt), 
S n (t) = sm(y/2(2n-l)gt), 



and the symmetric state: 



\4> + )ab= (|0) a |l) 6 +|l) o |0) 6 )/V2. 

So, the <£(p a b) concurrence of the (j2J density operator is given by 



when it is positive and is zero otherwise. Both terms on the right side hand of Eq. J3Jl are zero for n = whereas 
the second term is also zero for n = 1 and, for other values of n, the second term always reduce the concurrence. 
Therefore, as a function of the number n of excitations, the concurrence © acquires its maximum value for n = 1 at 
any time instant and it is given by €(p a b) — sin 2 (\/2.g£). On the other hand, for n ^> 1 the concurrence behaves 
as 



(4) 



3 



.... J 

2 4 6 8 10 

gt 

FIG. 1: Evolution (gt) of the concurrence for initially unexcited a and b atoms and the field in an Fock state. Black means 
value 1, maximum entanglement, white means value zero, and greys mean partial values of entanglement. 

The behavior of €(p a b) ~ 1 in the limit n ^S> 1 was found numerically by Tessier et al. [2^]. Figure ^shows the 
concurrence as a function of the n initial Fock state, and the gt adimensional time. Black means value 1, maximum 
entanglement, white means value zero, whereas greys mean partial values of entanglement. It can be seen that 
maximun value 1 is only reached for the initial condition |n = 1). In the next section we study how the concurrence 
is affected by the presence of an effective atomic environment when only one excitation is involved. 

III. EFFECTIVE HAMILTONIAN DESCRIPTION 

We consider a collection of A + 2 non-identical two-level atoms interacting with a single mode of a quantized field 
in an ideal cavity. Two atoms, labelled by subindexes a and 6, are resonant with the mode whereas the other A atoms 
interact dispersively with the mode. The Hamiltonian which drives the unitary dynamics of the whole system under 
the rotating wave approximation has the form 

A 

i—a,b J — 1 

A 

+ Sjffls+j+aU-j), (5) 

where and a are the usual one mode field operators, and s Zt ±j are the z components of the Pauli operators 
corresponding to the j— th two-level atom (j = a, b, 1, . . . , A). Atomic operators obey the standard SU(2) commutation 
relations, [s+,-, S-j] — 2s Z iSij and [s Z j,s±j] = ±s±j#y. Since the total number of excitations, represented by the 
operator N = a^a + y^_ n b 1 s Z j, is an integral of motion, the above Hamiltonian can be rewritten as follows: 

H = Lj f N + H int , (6) 

with 

A 

H m t = AjSzj + ^ gi{as +l + ata-i) 

j — 1 i—a : b 
A 

+ X^'( as +^' + atfi -i)' ( 7 ) 

3=1 

where — ujk — ujf, k = 1, . . . , A are the detunings between the transition of the &th atom and the mode frequency. 
Now, we assume that all A atoms are far from the resonance, so that Aj ^> j — 1,2,..., A. The effective 
Hamiltonian, approximately describing the interaction process, can be obtained from the interaction Hamiltonian 
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Q by using the method of Lie rotations |28l l29j , namely by applying to the Hamiltonian JSJ the following unitary 
transformation 

A 

V = e z, B = J2e j (as +j -ah- j ), (8) 
i=i 

where ej = gj/Aj <C 1. Neglecting terms of order higher than e 2 , we obtain the following effective Hamiltonian: 



A A 

He ff = y^( A j + fl[3 e 3'(l + 2a t a))s 2J - + A 9t(as +i + a f s_ 4 ) + - ^ j. 



j'=i 



2 

3>i=l 



(9) 



where we have defined the operator A = 1 + Sj=i e "j s zj- 

The last term in JjJ} represents an effective dipolar interaction between the non-resonant atoms, and its contribution 
to the system dynamics strongly depends on the internal resonance condition between atomic frequencies. Let us 
consider randomly distributed frequencies, such that they satisfy the condition Aj — Aj 3> g, i,j = 1,2, ..A. Then, 
the terms i ^ j in the last sum of the effective Hamiltonian © rapidly oscillate and can be neglected. Finally, the 
effective Hamiltonian, up to a constant energy shift, becomes: 

A A 

H eff = Y A 3 Sz 3 + ( X + 2fltfl ) Y 9i e 3 S *j 
3=1 3 = 1 

+X ^ 9i(as +i +a t s_ i ). 

i— a,b 

In the given approximation the total number of excitations "stored" in the non-resonant atoms is a constant of motion, 
which reflects a dispersive character of interaction. The first term in the above equation represents just transition 
frequency shifts of the non-resonant atoms, and commutes with the rest of the terms (so, it can be taken out of the 
Hamiltonian). The second term is the dynamic Stark shift and its contribution to the resonant dynamics, described 
by the last term, strongly depends on the state of the non-resonant atoms, which can be considered as a kind of 
atomic environment. 

Since the maximum entanglement in the system of two resonant atoms is reached when the total number of excitation 
is one, N = 1, we consider exclusively this situation. So, under the constraint that there is only one photon, the 
corresponding evolution operator can be found and, in the standard tensor product basis, it is given by 



U(t) = 



e lSt / 2 i„ + i -igXaLne 1 ^/ 2 -igXah^l 2 

-ig%^L n+x e m,i % y n - 1 

-igX^L n+x e^l 2 Y n - 1 Y n 

-igXatLn+xe-W 2 -igXa^L^e-^l 2 





-igXaLne-^' 2 
-igXaLne-^' 2 

e -m/2Az 



(10) 



where we have defined the operators: 



Y n 

An 



sin £l n t 



i(e^/ 2 i; l + e -^/ 2 i„ +1 ) 



cos Vt n t + i — L r , 



and the Rabi frequencies (l n depend on the field and on the environment variables as follows: 



fin 



+ 2g 2 X 2 n 



1/2 



A 

Y93 e 3 S z 
3 = 1 



So, the dynamics depends on the distribution of the different Rabi frequencies which appear as a contribution of A 
non-resonance distinguishable two-level atoms and it also depends on the initial state of the whole system. 
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FIG. 2: Evolution of the concurrence for initially unexcited a and b atoms and the field in the one-photon Fock state, with 
A = 7, g/g = 1, and g/A = 0.1. 

IV. EVOLUTION OF THE ENTANGLEMENT 

Now, let us suppose that the environment atoms are prepared in a coherent superposition of excited and ground 
states: 



I* 



0/ env 



I[j = i^M£_t!!k_L_£ £ . . .j> + i/ 2)l | S2 + 1/2 )2 . . . \ SA + 1/2) A , 



(ii) 



where the last sum is taken over all possible binary vectors s — {s±, S2, . . . sa}, Sj — s Z j = ±1/2. We capture the key 
features of the concurrence evolution by considering two particular cases of initial conditions. First, we suppose that 
the resonance atoms are initially in the ground state and that the field is in the one photon Fock state: 



l*> = |0).|0) 6 |l)/|*o) 



(12) 



Applying the evolution operator (|10fl to the i|12f) state and tracing up over the field and the off-resonance atomic 
environment, we obtain for the resonant atoms the following reduced density operator: 



Pab 



2 

^{sin 2 (^ m) t)|^ + )a fc a fc (^ + |+cos 2 (r!( m) i)|0) Q |0) fcQ (0| fc <0|}, 



2 A 



(13) 



with 



= |(^ m) ) 2 = |(E^5 m))2 ' 



A 



(m) 



(14) 
(15) 



where s^ m ^> is a certain arrangement of the s vector and Sj = ±1/2. We have neglected order corrections higher 
than or equal to e on the amplitudes. From now on the coupling constants for field-environment atoms are taken to 
be equal, gj = g for all j. The concurrence corresponding to the density matrix l|13|) takes the form 



2 

^) = ^Esin 2 (^t). 



(16) 



If the resonant atoms are initially prepared in the symmetric one excitation state and the field is in the vacuum 
Fock state, then the initial state of the whole system is the following tensor product: 



W= I°).H)>4_|1).I°), |0)/W< 



(17) 
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In a similar way as described in the previous case 1|12[) the concurrence takes the form 



2- 

£M = ^lE co S 2 (n[ m) t). (18) 



2 A 

m—l 

It is worth noting that the concurrences I|l(jll8|) are composed of many Rabi frequencies, which leads to a structure 
similar to the collapses and revivals in the Jaynes-Cummings model OH^- Nevertheless, in the present case the 
set of different frequencies is due to the presence of the dispersive atomic environment in contrast to the standard 
JCM where different Rabi frequencies appear as contributions of different Fock field states (recall that a well-defined 
collapse-revival structure requires a significant number of excitations |32l l33| ) . 

Let us consider randomly distributed numbers Aj, j = 1, . . . ,A with the mean A and the standard deviation a a- 
Then, the numbers e^, e 2 have the following mean values and standard deviations: 

e *A{ 1+ fr)> ''"W (19) 
~2 ~ 5 2 (. , „ 2g 2 a A .... 

e " * aH + ^ 

First, we will find the distribution of the l|14|l quantities. We assume that er e <§; e and a e i <C e 2 . Then, there are A + 1 
peaks corresponding to different values of the number of positive components of s, k = 0, . . . A; for a given value of 

fc there are C A — A\/[k\(A — fc)!] values of (e • s) and (e 2 ■ s) which are normally distributed in accordance with the 
central limit theorem. For the kth peak, the mean value and the standard deviation are given by 



k(A- 


-fc) 




1 ' 




-fc) 




- 1 



((?• s)) k = (k-A/2)e, a e , k = a. 

((?•§))* = (fc-A/2)?, * e 2, k =a £ 2 

Note that the first and the last peaks are infinitely narrow. 

The Rabi frequency distribution has A + 1 peaks (now the summation in (|16l I18J) is from fc = to A) , and the 
frequency corresponding to the fcth peak can be approximated as follows: 

^^V2 9 [i + ^(k-^)(i + ^- 2 (k-^))]. 



Thus, the expressions Ijltjfl and (|18|l for the concurrence can be approximated as follows: 

A 

2 A 



1 A 



fc=0 



(21) 



The d k separation between the fcth and the (fc + l)th peaks, and the 6 k width of the fcth peak are 

d k = v^[£ + 5 J^(2fc-A + l)], 

S k = V2g[l + (fc-|) + |^(fc-^)> e2 , fc . 

Then, considering the approximation of narrow peaks, S k < 4, i.e. cta *C A/v^4, the sum in (|21|l can be represented 
as a sum of Gaussians, that is 



with the width 



£(p a6 )«-± ^^Re > e ^ , (22) 

yH ' 2 2~g 2 ^2A<j ^ 

a fc — — CO 



, A 4 V2A 2 
a= WA- U ^ (23) 
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which grows with time. The 1221) sum reveals the collapse-revival structure of the concurrences l|16|) and (|18fl . The 
first collapse happens when 



and it is followed by revival at time: 



A 2 

gt R ~nk-^, k = l,2,.... (25) 

In Figure|21we show the exact evolution of the concurrence for initially unexcited atoms and the field in the one-photon 
Fock state in the presence of the environment atoms. One can observe that the entanglement also reaches its maximum 
value. We can estimate from (|24|) and (|25|l the time scale required to observe the environment induced collapse-revival 
structure. Taking the typical values of the interaction constant from (25|: g/2-K — 2AkHz and K/2ir ~ 70kHz, we 
obtain i# ~ 100/is which is of order of the passage time of the atom through the cavity (cold atoms, v ~ lOOm/s) 
and less than the photon lifetime ~ 160/is. The collapse time is y/A times less than tn. 



V. CONCLUSIONS 



In summary, we have studied the dynamics of the concurrence of two atoms resonantly interacting with a cavity 
mode in the presence of many off-resonance atoms. We have shown that, for random distribution of atomic detunings 
and initially symmetric excitation of non-resonant atoms, the coherent influence of the environment can be separated 
from the dephasing. The coherence influence of the environment reflects in the appearance of the collapse-revival 
structure of the concurrence, with an average value of one half. This behavior is induced only by the presence of the 
dispersive environment. It is worth noting that the entanglement of formation in revival periods has as extreme value 
1, which means that the state of the bipartite system becomes pure at those times. 
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